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We study the N/Z and N/A dependence of balance energy i^hai) for isotopic series 
of Ca having N/Z (N/A) varying from 1.0 to 2.0 (0.5 to 0.67). We show that the N/Z 
(N/A) dependence of "Ebai is sensitive to symmetry energy and its density dependence at 
densities higher than saturation density and is insensitive towards the isospin dependence 
of nucleon-nucleon (nn) cross section and Coulomb repulsion. We also study the effect 
of momentum dependent interactions (MDI) on the N/Z (N/A) dependence of 'Ej^ai- We 
find that although MDI influences the Yj^ai drastically, the N/Z (N/A) dependence of Efea; 
remains unchanged on inclusion of MDI. 
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1 Introduction 



With the availabiUty of high-intensity radioactive beams at many facihties as well as a 
number of next generation beam facilities being constructed or being planned [1,2], the 
studies on the role of isospin degree of freedom have recently attracted a lot of attention 
in both nuclear physics and astrophysics. The ultimate goal of such studies is to extract 
information on the isospin dependence of in-medium nuclear effective interactions as well 
as equation of state (EOS) of isospin asymmetric nuclear matter. The later quantity 
especially the symmetry energy term is important not only to nuclear physics community 
as it sheds light on the structure of radioactive nuclei, reaction dynamics induced by rare 
isotopes but also to astrophysics community as it acts as a probe for understanding the 
evolution of massive stars and supernova explosion [3]. It is worth mentioning that the 
equation of state of symmetric nuclear matter has been constrained up to densities 5 times 
the normal nuclear matter density through the measurements of transverse flow as well as 
its disappearance along with other collective flows (like radial flow, elliptic flow) [4] and 
of subthreshold kaon production in relativistic nucleus-nucleus collisions [5] . 

Although the nuclear symmetry energy at normal nuclear matter density is known 
to be around 30 MeV [6], its values at other densities are poorly known. Heavy-ion 
collisions induced by radioactive beams provide unique opportunities to investigate the 
isospin-dependent properties of asymmetric nuclear matter, particularly the density de- 
pendence of symmetry energy [7]. Experimentally symmetry energy is not a directly 
measurable quantity and has to be extracted from observables related to symmetry en- 
ergy. Over the last decade a large number of experimental observables have been proposed 
like neutron/proton ratio of emitted nucleons [8], the neutron-proton differnetial flow [9], 
the t/^He [10], 7r~/7r+ [6,11], E-/E+ [12], and K7ii'+ [13] ratios and so on. A recent 
analysis of data has led to a symmetry energy term of the form Egym — 31-6(^)^ MeV 
with 7 = 0.4-1.05 for densities between O.lpo and 1.2po [14]. However, for all the above 
mentioned observables the Coulomb force of charged particles plays an important role. 
It competes strongly with symmetry energy. Recently Gautam and Sood [15] studied 
the relative contribution of Coulomb force and symmetry energy in isospin effects on the 
collective transverse flow as well as its disappearance for isobaric systems throughout the 
mass range and colliding geometry. They clearly demonstrated the dominance of Coulomb 
repulsion over the symmetry energy. The collective transverse in-plane flow disappears at 
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a particular energy called as balance energy [15, 16]. In recent communication, Gautam 
et al. [17] has studied the transverse momentum for a neutron rich system ^^Ca+^^Ca in 
the Fermi energy as well as at high energies. There they find that transverse momentum 
is sensitive to the symmetry energy as well as its density dependence in the Fermi energy 
region. Motivated by those results we here study the E^ai as a function of N/Z and N/A 
of the system for an isotopic series. We here choose the isotopes so that the Coulomb 
repulsion is same for all the systems, since as mentioned previously that Coulomb plays 
much dominant role as compared to symmetry energy in isospin effects. Here we will 
demonstrate that the N/Z (N/A) dependence of Em for the isotopes of same element is 
a sensitive probe for the symmetry energy as well as its density dependence. To check 
the sensitivity of N/Z (N/A) dependence of towards density dependence of symmetry 
energy, we have calculated the Ebai throughout the isotopic series for different forms of 
symmetry energy Fi{u), F2{u), and F^lu), where u — The different forms are de- 
scribed later. The present study is carried out within the framework of Isospin-dependent 
Quantum Molecular Dynamics (IQMD)[18] Model. Section 2 describes the model in brief. 
Section 3 explains the results and discussion and section 4 summarizes the results. 

2 The model 

The IQMD model treats different charge states of nucleons, deltas and pions explicitly, as 
inherited from the Vlasov-Uehling-Uhlenbeck (VUU) model. The IQMD model has been 
used successfully for the analysis of a large number of observables from low to relativistic 

energies. One of its versions QMD model has been quite successful in explaining various 
phenomena such as multifragmentation [19], collective flow [20], and hot and dense nuclear 
matter [21] as well as particle production [22]. The isospin degree of freedom enters into 
the calculations via symmetry potential, cross sections and Coulomb interaction. 
In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

fi{r,P,t) = ^exp(-[r-rl(t)]2— ) x exp(-[p - pl(t)]2— ) (1) 

Nucleons are initialized in a sphere with radius R = 1.12 A^/"^ fm, in accordance with 
liquid-drop model. Each nucleon occupies a volume of h^, so that phase space is uni- 
formly filled. The initial momenta are randomly chosen between and Fermi momentum 
(Pf)- The nucleons of the target and projectile interact by two- and three-body Skyrme 
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forces, Yukawa potential, Coulomb interactions, and momentum-dependent interactions 
(MDI). In addition to the use of explicit charge states of all baryons and mesons, a symme- 
try potential between protons and neutrons corresponding to the Bethe-Weizsacker mass 
formula has been included. The hadrons propagate using Hamilton equations of motion: 

dfi d{H) dpi d{H) 



dt dpi ' dt dfi 



(2) 



with 



(H) = {T) + {V) 

i ^ i j>i 

xfj(r',p',t)dfdf' dpdp'. (3) 



The baryon potential , in the above relation, reads as 



V' '/ ' Skyrme ' Yukawa ' Coul ' mdi ' * sym 

= M(r'-r)+M(f'-rV-^(^^)] 
,^ exp(|(r^-f)|//x) , Z,Z,e^ 

+^-3 — TTT^i -^M / X r 



(|(f'-f)|//x) \{f'-r)\ 
+tihi^[t^{p' -p)^ + l]5{r' -f) 

+te-TsiTs,d{ri' -fj). (4) 

Here Zi and Zj denote the charges of ith and jth baryon, and T^i and T^j are their 
respective components (i.e., 1/2 for protons and —1/2 for neutrons). The parameters n 
and ti,....,tQ are adjusted to the real part of the nucleonic optical potential. For the density 
dependence of the nucleon optical potential, standard Skyrme-type parametrization is 
employed. We also use the isospin and energy- dependent cross section a = 0.8 (J^^^^. The 
details about the elastic and inelastic cross sections for proton-proton and proton-neutron 
collisions can be found in [18,23]. The cross sections for neutron-neutron collisions are 
assumed to be equal to the proton-proton cross sections. Explicit Pauli blocking is also 
included; i.e. Pauli blocking of the neutrons and protons is treated separately. We assume 
that each nucleon occupies a sphere in coordinate and momentum space. This trick yields 
the same Pauh blocking ratio as an exact calculation of the overlap of the Gaussians 
will yield. We calculate the fractions Pi and P2 of final phase space for each of the two 
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scattering partners that are already occupied by other nucleons with the same isospin as 
that of scattered ones. The colhsion is blocked with the probability 

Pbiock = 1 - [1 - mm(Pi, 1)][1 - mm(P2, 1)], (5) 

and, correspondingly is allowed with the probability 1 - P block- For a nucleus in its 
ground state, we obtain an averaged blocking probability (Puock) — 0.96. Whenever an 
attempted collision is blocked, the scattering partners maintain the original momenta 
prior to scattering. The different forms of symmetry energy are obtained by changing the 
density dependence of the potential part of the symmetry energy (last term in eq. (4)). 
The various forms are Fi{u) oc u, F2{u) oc vP^^, F^^u) oc (where u = -^). F4 represents 
calculations without symmetry potential. 

3 Results and discussion 

We have simulated several thousand events at incident energies around balance energy 
in small steps of 10 McV/nucleon for each isotopic system of Ca+Ca having N/Z (N/A) 
varying from 1.0 to 2.0 (0.5-0.67). i.e. Ca^o+Ca^o, Ca^^+Ca^^ Ca^8+Ca^^ Ca^^^Ca^^ 
Ca^^+Ca^^, and Ca®°+Ca^° for the semicentral colliding geometry range of 0.2 - 0.4. Such 
systematic studies performed at low incident energies using various fusion models have 
shown a linear enhancement in the fusion probabilities with neutron content [24]. We 
use soft equation of state with and without MDI, labeled respectively as Soft and SMD. 
The calculations with this choice of equation of state and cross section were in good 
agreement with the data throughout the coUiding geometry [25]. The IQMD model has 
also been able to reproduce the other data (example, high energy proton spectra, gamma 
production) at incident energies relevant in this paper [26, 27]. The reactions are followed 
till the transverse flow saturates. The saturation time is around 100 fm/c for the systems 
in the present study. For the transverse flow, we use the quantity "directed transverse 
momentum {pf^)" which is defined as [28, 29] 

1 ^ 

(Pf ) = jJ2^'3n{yii)}Px{i), (6) 

i=l 

where y{i) is the rapidity and Px{i) is the momentum of i^^ particle. The rapidity is 
defined as 

y(.) = lln^W + Pf;. (7) 
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Figure 1: (Color online) Em as a function of N/Z (upper panel) and N/A (lower panel) of 
system for Egym oc Fi{u) and F4. Lines are linear fit proportional to m. Various symbols 
are explained in the text. 



where E(z) and P2(z) are, respectively, the energy and longitudinal momentum of i*^ 
particle. In this definition, all the rapidity bins are taken into account. A straight line 
interpolation is used to calculate the Ei^ai- It is worth mentioning that the Ei^ai has the 
same value for all fragments types [31-34]. 

In fig. 1(a) we display the E{,ai as a function of N/Z of the system. Solid green 
circles represent the calculated Ehaj. Lines are the linear fit to E^a;. We see that E;,a; 
follows a linear behavior oc m*N/Z. As the N/Z of the system increases, the mass of the 
system increases due to addition of neutron content. In addition, the effect of symmetry 
energy also increases with increase in N/Z. To check the relative contribution of increase 
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in mass with N/Z and symmetry energy towards the N/Z dependence of E^ai, we make 
the strength of symmetry energy zero and calculate Ebai- The results are displayed by 
open circles in fig. 1(a). E^ai again follows a linear behavior oc m*N/Z. However, Ef,ai 
decreases very slightly with increase in N/Z, whereas when we include the symmetry 
energy also in our calculations then the | m \ increases by 3 times which shows that N/Z 
dependence of E^ai is highly sensitive to the symmetry energy. The shght decrease in the 
Ej,ai with N/Z (for calculations without symmetry energy) is due to increase in number 
of nucleon-nucleon collisions. To further explore this point, we switch off the symmetry 
energy and also make the cross section isospin independent (i.e. (T„p = (T„„ and calculate 
Ebai for two extreme N/Z. The results are displayed in fig. 1(a) by open squares. Again 
Ebai follows a linear behavior. We see that the E^a; for both ^°Ca + *^°Ca and ^°Ca + 
^°Ca increases as expected. However, the increase in Ebai for system with N/Z = 1 is 
more as compared to the system with N/Z = 2. This is because with increase in N/Z 
the neutron number increases due to which neutron-neutron and neutron-proton collisions 
pairs increase. However, the increase in number of neutron-neutron collision pairs is much 
larger as compared to neutron-proton collision pairs. Therefore, the possibility of neutron- 
proton collision is much less in system with N/Z = 2. That is why the effect of isospin 
dependence of cross section decreases with increase in N/Z. 

In fig. 1(b), we display the Ebai as a function of N/A of the system. Symbols have 
same meaning as in fig. 1(a). Again Ebai follows a linear behaviour with m = -191 and 
-68, respectively, for Fi (u) and F4. However, the percentage difference AEbai % (where 
AEbai % — ''"'f^ („)'"'' ) is same (about 65%) in both the figs. 1(a) and 1(b) which shows 
that the effect of symmetry energy is same whether we discuss in terms of N/Z or N/A. 

As stated in literature, the isospin dependence of collective flow and its disappearance 
has been explained as the competition among various reaction mechanisms, such as, nn 
collisions, symmetry energy, surface property of the colliding nuclei and Coulomb force 
[30]. Here we aim to show that the N/Z and N/A dependence of Ebai is sensitive to the 
symmetry energy only. Since we are taking the isotopic series, the effect of Coulomb will 
be same for all the reactions. We have also checked that the N/Z dependence of Ebai is 
insensitive to the EOS of symmetric nuclear matter. Moreover, as mentioned previously, 
the equation of state of symmetric nuclear matter has been constrained up to densities 
five times the normal matter density. In the present case as the N/Z of the system 
increases, the number of neutrons also increases. Since we are using isospin-dependent 
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Figure 2: (Color online) E^a/ as a function of N/Z (left panels) and N/A (right panels) of 
system for isospin independent cross section for (a) and (c) N/Z = 1.0 to 2.0 (b) and (d) 
N/Z = 1.2 to 2.0. Lines are linear fit proportional to m. Various symbols are explained 
in the text. 



nn cross section, so to check the sensitivity of N/Z and N/A dependence of E^a/ to the 
isospin dependence of cross section, we calculate the Efea^ throughout the isotopic series 
by making the cross section isospin independent (fig. 2(a) open orange triangles, left 
panels). Again E^ai follows a linear behavior with m = -40. We find that although E^a; for 
individual system is very sensitive to the isospin dependence of cross section. However, 
N/Z dependence of Efe^; (for isotopic series) is much less sensitive to the isospin dependence 
of cross section. 

In fig. 2(b) (left panels), we show E^ai as a function of N/Z of the system for the N/Z 
range from 1.2 to 2.0. We find that the sensitivity of N/Z dependence of E^a; towards the 
isospin dependence of cross section decreases further. Now m is -34 (-37) for calculations 
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Figure 3: (Color online) Em as a function of N/Z (upper panel) and N/A (lower panel) 
of system for Egym. oc F2{u) and Fs{u). Lines are linear fit proportional to m. Various 
symbols are explained in the text. 



with (without) isospin dependence of cross section. Thus the N/Z dependence of for 
neutron-rich isotopes is sensitive only to symmetry energy. In figs. 2(c) and 2(d) (right 
panels) we display similar plots as in the corresponding left panels but now we plot E^ai as 
a function of N/A of the system. Again the percentage difference between the two curves 
in both upper panels is same (about 18 %) and same in lower panels as well (about 8 %). 

In fig. 3 (a) (3b) we display the N/Z (N/A) dependence of E^ai for different forms 
of symmetry energy; Fi(u) (solid circles), F2(u) (diamonds), and F3(u) (pentagons). For 
all the cases E^ai follows a linear behavior. Clearly, N/Z (N/A) dependence of Ehai is 
sensitive to the density dependence of symmetry energy as well. For a fixed N/Z (N/A) 
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Figure 4: (Color online) The time evolution of < p^*** > for different forms of symmetry 
energy for different bins at b/bmax=0.2-0.4 . Lines are explained in the text. 



stiff symmetry energy Fi(u) shows less E^ai as compared to soft F2(u) whereas super stiff 
symmetry energy F3(u) shows more E^ai as compared to F2(u). 

To explain the above mentioned feature, we calculate for ^''Ca + ^°Ca the transverse 
flow of particles having p/po < I (denoted as BIN 1) and particles with p/po > 1 (denoted 
as BIN 2), separately at all time steps for symmetry energy Fi(u), F2(u), and F3(u). The 
incident energy is taken to be 100 MeV/nucleon. The results are displayed in fig. 4. Solid 
(dashed) lines represent the pf*^ of particles lying in BIN 1 (BIN 2). Dotted line represent 
the total < pf^ >. We see that the total < pf^ > is maximum for stiff symmetry energy 
and minimum for super stiff symmetry energy. During the initial stages of the reaction, 
< pf^ > due to particles lying in BIN 1 remains positive for Fi(u) and F2(u) because in 
the spectator region, repulsive symmetry energy will accelerate the particles away from the 
overlap zone. The effect is more pronounced for Fi(u) as compared to F2(u). Moreover, 
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for Fi(u) and F2(u), this interval is about 5-25 fm/c and 5-20 fm/c, respectively. This 
is because although for F2(u), the effective strength of symmetry energy will be more 
for low density particles as compared to Fi(u), however, in the central dense zone the 
effective strength of F2(u) will be less i.e. in the central dense zone, F2(u) will be less 
repulsive, therefore for F2(u), there will be more attractive force on the particles lying in 
the spectator region towards the central dense zone as compared to that in case of Fi(u). 
That is why during the initial stages the peak value of < pf^ > as well as the duration 
for which it remains positive is less for F2(u) as compared to Fi(u) (compare shaded area 
in fig. 4(a) and 4(b). This decides the value of < p^*'' > at saturation, which is more 
for Fi(u) as compared to F2(u). In case of F3(u) (fig. 4(c)) for particles lying in BIN 
1, i.e. (p/po ^ 1 ); the strength of symmetry energy will be much smaller which is not 
sufficient to push the particles away from the overlap zone. Therefore, the < pf"' > of 
BIN 1 particles remains zero during the initial stages. This leads to least value of final 
state < p'^'^ > for super stiff symmetry energy as compared to stiff and soft symmetry 
energy. The < p*^ > due to particles in BIN 2 (dashed line) decreases in a very similar 
manner for all the different symmetry energies between 0-10 fm/c. Between 10-25 fm/c, 
< p"^^ > for ¥^{u) decreases more sharply as compared to in case of Fi(u) and F2(ii). This 
is because in this time interval the particles from BIN 1 enters into BIN 2 and < p'^^ > 
of particles entering BIN 2 from BIN 1 in case of Fi(-u) and F2{u) will be less negative 
due to stronger repulsive symmetry energy as compared to in case of F3 (u) (see Ref. [17] 
also). During the expansion phase, i.e. after 30 fm/c the total < pf^ > and < pf^ > of 
BIN 1 particles overlap as expected. Therefore, the effect of symmetry energy on the low 
density particles during the initial stages decide the fate of the final < pf'^ > and hence 

Since one cannot use radioactive isotopes as targets, therefore, as a next step we fix 
the target as a stable isotope ^°Ca and vary the projectile from ^°Ca to ^°Ca and calculate 
E''"'. In this case the N/Z (N/A) of the reaction varies between 1 to 1.5 (0.5 to 0.6) and 
the asymmetry 5 = of the reaction varies from to 0.2. The results are displayed 

by solid green stars in figs. 5(a) and (b) (upper panels). The solid green circles represent 
the calculations for symmetric reactions with N/Z (N/A) varying from 1 to 2 (0.5 to 
0.67), i.e. '^°Ca+'^°Ca to ®°Ca-|-^°Ca. Lines represent the linear fit oc m. We see that 
N/Z (N/A) dependence of E^a; is same for both the cases. We also find that when we use 
stable target ^°Ca and radioactive target ^°Ca, the N/Z (N/A) decreases from 2 (0.67) 
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Figure 5: (a) and (b) (upper panels) E^ai as a function of N/Z (left panel) and N/A (right 
panel) of system for E^y^ oc Fi{u) with ''^Ca as target (stars), (c) and (d) (lower panels) 
Ebai as a function of N/Z (left panel) and N/A (right panel) of system for Egym Fi{u) 
with SMD EOS (left triangles). Circles represent the values of Ef,ai as in fig. 1. Lines are 
linear fit proportional to m. 



in case of ^°Ca+®°Ca to 1.5 for ^"CaH- Ca, so the Ehai also decreases. Now the Ehai for 
^'^Ca+'^°Ca has same value as in case of symmetric reactions with N/Z (N/A) = 1.5 (0.6) 
i.e. the value of Ebai is decided by the N/Z (N/A) of the system and is independent of 
the asymmetry of the reaction in agreement with [35]. 

It has also been reported in literature that the MDI affects drastically the collective 
flow as well as its disappearance [36]. To check the influence of MDI on the N/Z (N/A) 
dependence of Ebai we calculate the Ebai for the whole N/Z (N/A) range from 1 to 2 (0.5 
to 0.67) for the symmetric reactions with SMD equation of state and symmetry potential 
Fi{u). The results are shown in figs. 5(c) and (d) (lower panels) by solid left triangles. 
We find that although the MDI changes drastically the absolute value of E^ai ( by about 
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30%), however the N/Z (N/A) dependence of E^ai remains unchanged on inclusion of 
MDI. Therefore, the dependence of E^ai as a function of N/Z (N/A) on the symmetry 
energies of other different forms (F2(m) and F3(m)) is also expected to be preserved on 
inclusion of MDI. 

4 Summary 

We have shown that the N/Z (N/A) dependence of Ef^ai for the isotopic series of Ca+Ca 
is a sensitive probe to the symmetry energy as well as its density dependence at densities 
higher than saturation density and is insensitive to other isospin effects like Coulomb 
repulsion, and isospin dependence of nucleon-nucleon cross section. We have also studied 
the effect of MDI on the N/Z (N/A) dependence of E^ai. We find that although MDI 
influences the Ehai drastically, the N/Z (N/A) dependence of Ei^ai remains unchanged on 
inclusion of MDI. 

This work has been supported by a grant from Indo-French Centre For The Promotion 
Of Advanced Research (IFCPAR) under project no. 4104-1. 
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